SiMilirS DDtI original research article 

published: 05 November 2013 
doi: 10.3389/fncom. 2013.00145 



COMPUTATIONAL NEUROSCIENCE ^ ^ °' '° ^ ° ' 




A possible role of astrocytes in contextual nriemory 
retrieval: An analysis obtained using a quantitative 
framework 

Shivendra Tewari^-^* and Vladimir Parpura^** 

' Biotechnology and Bioengineehng Center, Medical College of Wisconsin, Milwaukee, Wl, USA 
' Department of Physiology, Medical College of Wisconsin, Milwaukee, Wl, USA 
^ Department of Neurobiology, University of Alabama at Birmingham, Birmingham, AL, USA 
* Department of Biotechnology, University of Hijeka, Rijeka, Croatia 



The hippocampus is central to our understarnding of memory formation and retrieval. 
Its CA1 region is known for encoding contextual memory. Here, using a computational 
approach, which embeds existing physiological data, we propose a particular role of 
astrocytes in contextual memory retrieval. We provide a quantitative framework under 
which the astrocyte modulates the firing of a context-associated CA1 pyramidal neurons, 
resulting in a prominent tuning of neurons to a delta rhythm. Using the very framework, we 
further studied astrocytic function in the modulation of neuronal firing under pathological 
conditions, i.e., during astrocytic induction of epileptiform discharge in CA1 pyramidal 
neurons. Thus, we provide a quantitative framework that would aid understanding of the 
Schaffer collateral-CAl tripartite synapse in health and disease. 
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INTRODUCTION 

One of the most intriguing properties of the brain is to learn and 
form new memories. Memories are essential for guiding effective 
actions, for understanding the words we read, for recognizing the 
objects we see, etc. It is well established that retrieval of a particu- 
lar memory is facilitated when the context at the time of memory 
retrieval is similar to that of the initial learning experience (Maren 
and Holt, 2000). This phenomenon, termed as context-dependent 
memory, also provides an explanation as to why memory recall 
is impaired when people, place or things are experienced "out 
of context." The process by which the context is embedded into 
learning of information, and can later facilitate memory retrieval 
is known as contextual memory retrieval (Maren and Holt, 2000). 
Pavlovian conditioning, i.e., classical conditioning, is an exemplar 
form of contextual fear memory learning and retrieval (Maren, 
2001). 

Ever since memory deficits, initially characterized in patients 
PB and FC (Milner and Penfield, 1955) and later HM (Scoville 
and Milner, 1957), were associated with hippocampus damage, 
it has become clear that the hippocampal region plays a crucial 
role in the acquisition of new semantic and episodic memory 
(Andersen et al., 2007); but see (Corkin, 2002). Apart from stor- 
ing memories, the appropriate retrieval of stored memories is also 
of paramount importance for maintaining relationships, achiev- 
ing personal and professional goals, and carrying out even simple 
day to day activities. The critical role of the hippocampus in 



memory retrieval has gone largely uncontested and has been 
supported by various studies including lesion studies (Holt and 
Maren, 1999), unit studies (Paller and McCarthy, 2002), as well 
as functional neuroimaging (Eldridge et al., 2000). 

Lee and Kesner (2004) reported a sub-region specific contri- 
bution of the hippocampus in memory acquisition and retrieval. 
Interpretation of the result suggests that the CAl region is 
uniquely associated with a reliable retrieval of the contextual 
memory. Such a role for the CAl region in contextual memory 
retrieval has also been highlighted in several studies (Hall et al., 
2001; Khuchua et al, 2003; Murchison et al, 201 1). Previous work 
suggests that the firing synchrony in neighboring CAl pyramidal 
neurons might be a putative mechanism of contextual memory 
retrieval in the CAl region (Hall et al., 2001; Siekmeier et al., 
2007; Takahashi and Sakurai, 2009). However, data are accu- 
mulating suggesting that the coordination between neurons and 
astrocytes may play a role in information processing and memory 
retrieval. 

Fellin et al. (2004) demonstrated that the activation of Ca^+ 
elevation in astrocytes, by various stimuli, produces slow inward 
currents (SICs) in 2-12 CAl pyramidal neurons, synchronously. 
They found that repetitive stimulation of astrocytes leads to syn- 
chronous firing of the same CAl pyramidal neurons, mainly due 
to the activation of neuronal NR1/NR2B subunits of the extra- 
synaptic N-Methyl-D-aspartate receptors (eNMDARs). It should 
be noted that NMDAR NR2B subtype plays a crucial role in 
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FIGURE 1 I Schematic diagram of the proposed networit model of 
contextual memory retrieval within the hippocampus where an 
astrocyte modulates synchronous firing of a pack of four CA1 
pyramidal neurons during context-dependent memory retrieval. The 

pack of neurons is aligned within 100 |im of the stratum pyramidale along 
the input, i.e., Schaffer collateral (SC) axons of the CA3 pyramidal neuron, 
which makes synapses with the spines of the proximal apical dendrites of 
CAl pyramidal cells in the stratum radiatum. The field potential recording 
electrode is placed next to the CAl neuron 4, at the most distal end to the 
input of this pack. The fine processes of the astrocyte, allied with the 
neuronal pack and located in the stratum oriens, listen to the synaptic 
transmission which takes place in the stratum radiatum, at the synapse 
made by the CAS pyramidal neuron onto the lead CAl pyramidal neuron 1 
(the utmost left, proximal to the input). The astrocyte integrates received 
information and relays the signal to the following CAl pyramidal neurons 
within the pack (in this case the three remaining CAl pyramidal neurons 
2-4, distal to the input) associated with the context of the input signal from 
the CA3 pyramidal neuron concurrently, the astrocyte also potentiates the 
SC-CAl synapse strengthening, the small network activity. 



encoding and retrieval of contextual fear memory (Zhao et al., 
2005; Zhang et al, 2008; Zhuo, 2009); but see Nakazawa et al. 
(2004). Thus, it appears that the interaction/signaling between 
neurons and astrocytes plays a role in information processing and 
memory retrieval. 

Astrocytes in the hippocampus have numerous highly 
branched fine processes that extend to distances up to 100 [im and 
enwrap asymmetrical (mainly excitatory/glutamatergic) synapses 
of a relatively small group of neurons (Fellin et al., 2006); of 
note, individual astrocytic domains do not overlap (Bushong 
et al, 2002, 2004; Ogata and Kosaka, 2002; Nedergaard et al, 
2003; Oberheim et al, 2009). Thus, strategically, astrocytes 
are ideally positioned to sample electrical activity patterns of 
the Schaffer collateral input from CA3 pyramidal neurons and 
retrieve their associated context in neighboring CAl pyramidal 
neurons. During the retrieval, the firing pattern of the pyramidal 
neurons, in the CA3-CA1 region, is quite similar to that recorded 
during memory formation (Small et al, 2001). As astrocytes 
play an important role in short-term and long-term potentiation 
of individual synapses (Martin et al., 2007; Henneberger et al., 
2010), it is possible that they also contribute to various aspects of 
memory. We propose that astrocytes facilitate the memory storage 
and retrieval by sensing encoding/retrieval patterns of neurons 
that reside within individual astrocytic domains and, in turn, 
astrocytes signal the associated CAl pyramidal neurons. 

Based on these recent experimental observations, we devel- 
oped a computational framework/model, where an astrocyte 
modulates synchrony among CAl pyramidal neurons within a 
domain of 100 |xm in diameter, which qualitatively reproduces 
the biological experimental findings (Fellin et al., 2004). Power 
spectral density analysis of the computed local field potentials 
indicates a shift in neuronal firing rhythm from broad-frequency 
range in the absence of astrocytic gliotransmission to a delta- 
dominant rhythm when gliotransmission is turned on. As delta 
rhythm has been implicated in cognitive domain (Knyazev, 2012) 
and consolidation of declarative memories (Ruch et al, 2012), 
this might be a putative mechanism of memory retrieval tuned 
by astrocytes, which will require further experimental verification 
in biological preparations. Apart from a possible role in mem- 
ory retrieval, using the computational model, we propose that 
astrocytes may lead to epileptiform discharges in CAl pyrami- 
dal neurons, as has been previously demonstrated in biological 
preparations (Kang et al, 2005; Tian et al., 2005). Such patholog- 
ical discharges occuring under computational conditions warrant 
an increase in eNMDAR plasmalemmal activity that could arise 
from an overexpression of such receptors and/or positive mod- 
ulation of their conductance (Omkumar et al, 1996; Liao et al., 
2001). 

MATERIALS AND METHODS 
OVERVIEW OF THE MODEL ARCHITECTURE 

The complete network (shown in Figure 1) comprises of three 
major components: a CA3 pyramidal neuron, four CAl pyra- 
midal neurons and an astrocyte in the stratum oriens (s.o.). 
A CA3 pyramidal neuron is computationally described by Traub's 
branching dendrite model of the CA3 pyramidal neuron com- 
prising of 21 compartments (Traub et al., 1994). Each of the 



four CAl pyramidal neurons is described by the 19-compartment 
CAl pyramidal neuron model (Traub et al, 1991), which was 
modified to incorporate four additional compartments account- 
ing for the axon initial segment (IS), axon proper, and two 
spines (placed at the 6*'' and 12*'' compartment of the origi- 
nal CAl pyramidal neuron model) to yield a 23-compartment 
CAl pyramidal neuron model. Axon and spine compartments 
were added to form CAl-CAl synapses in the s.o. at the 
basal dendrites. The dendritic spine model utilized here is 
from Tewari and Majumdar (Tewari and Majumdar, 2012a), 
but was modified to include synaptic NMDARs and eNM- 
DARs. The astrocyte model is also adopted from Tewari and 
Majumdar (Tewari and Majumdar, 2012a,b), and it accounts 
for gliotransmitter (glutamate) release/astrocyte-dependent SICs 
through the activation of NR2B-containing receptors on the CAl 
neurons. It is assumed that the CA3 pyramidal neuron makes 
25 synapses with the CAl neuron 1 (Figure 1). The number 
of synapses at the CAl-CAl synapses is assumed to be 20. 
These numbers were obtained by a hit and trial approach, and 
reflect the smallest numbers possible to have a successful sig- 
nal transduction of the proposed network. All associated model 
equations and parameters are provided in the supplementary 
materials. 

As shown in Figure 1, the CA3 neuron forms synapses with 
the CAl neuron 1, at the spines in the stratum radiatum (s.r.). To 
model glutamate release from the CA3 and CAl neurons, we use a 
continuous function that transforms axon voltage into glutamate 
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concentration in the cleft (Destexhe et al., 1994): 
[g] 



1 + exp[- (Vaxon - Vp) /Kp]' 



(1) 



where ^max is the maximum concentration of glutamate in the 
synaptic cleft, Vaxon is the membrane potential in the axon 
compartment of the glutamate releasing neuron, Kp gives the 
steepness, and Vp is the axon voltage at which glutamate con- 
centration is half of the maximum (see Table 1). The glutamate 
released from CA3 neuron is sensed by the astrocyte in the s.o. 
(Perea and Araque, 2005) leading to astrocytic Ca^+ elevation. 
Simultaneously, glutamate released from the CA3 neuron leads to 
an inward current in the CAl neuron 1. The SICs evoked by the 
astrocyte through eNMDARs are: 



^SIC = gNRlB - h-Vs. 



(2) 



where ^nr2B is the conductance through NRl/NR2B-containing 
receptors, Vj is the spine membrane potential and h is the gating 
variable of eNMDAR that depends on the astrocyte as described 
by the following equation: 



dh 
dt 



a 



h 



(3) 



The above relation is regulated by astrocytic Ca^+ by a, which 
has a value of 0.4 (representing the maximum open probabil- 
ity of eNMDAR) whenever astrocytic Ca^+ goes above 200 nM 
(Parpura and Haydon, 2000; Nadkarni and Jung, 2003; Tewari 
and Majumdar, 2012b), otherwise set to zero. The time constants 
Ti and T2 govern the rise-time and decay- time of the eNMDAR 
currents, and are based on the experimental observation (Fellin 
et al, 2004); parameter values are listed in Table 1. 

A given CAl neuron forms synapses with the other CAl neu- 
rons in the basal dendrite region (Andersen et al, 2007); in the 
proposed network model, the synapses are formed at the 6* com- 
partment of the CAl neuron model. The sequence in which the 
CAl neurons make synapses is as follows: 1^2, 2— >-3, 3— >-4 
(Figure 1, left to right); the CAl neuron 4 is not connected to 
any other neuron. The astrocyte senses the signal (glutamate) 
released by the CA3 neuron and potentiates the above CAl neu- 
ron synapses by releasing gliotransmitter (glutamate) in the vicin- 
ity of the eNMDARs on the spines, which leads to synchronized 
firing of the CAl neurons. The astrocyte releases glutamate from 



Table 1 | List of parameters involved at the synapses. 



Symbol 


Value 


Reference 


9'max 


3.5-4 mM 


Danbolt, 2001 


Vp 


72 mV""" 


Destexhe et al., 1994 


Kp 


5mV 


Destexhe et al., 1994 


£'NR2B 


0.6 nS 


Gasparini et al., 2004 


ti 


92.3 ms 


Fellin et al., 2004 


12 


568.5 nns 


Fellin et al., 2004 



glutamate secretory vesicles utilizing Ca^+ -dependent exocytosis 
as soon as astrocytic Ca^+ concentration exceeds the thresh- 
old Ca^^ concentration for triggering secretion [see (Tewari and 
Majumdar, 2012a) for details]. 

COMPUTER SIMULATIONS 

The proposed network model invoked a total of 1247 ordinary 
differential equations which were solved using a quasi-constant 
step sized backward-difference numerical solver implemented in 
MATLAB (The Math Works, Natick, MA). The system of ordinary 
differential equations was solved for a duration of 50 s (model 
time) which took approximately 8 h (real time) to compute. 

LOCAL FIELD POTENTIAL CALCULATIONS 

The somata of the CAl neurons were spaced 25 \im apart. The 
recording electrode was placed 25 [im away from the soma of the 
CAl neuron 4 (see Figure 1). The local field potential was calcu- 
lated using line source approximation (Gold et al., 2006). Based 
on the ionic currents generated from the network simulation 
(in the presence and absence of the astrocyte), trans-membrane 
currents across each compartment were calculated for each of 
the CAl neurons. Please note that the contribution from IS and 
axon compartments was neglected because myelinated axon and 
nodes of Ranvier make no contribution to the local field poten- 
tial (Gold et al., 2006). Afterwards, linear source approximation 
was used to compute local potentials at the recording electrode 
using a custom-written MATLAB program. Mathematically, field 
potential is calculated using the following expression: 



3>(f) 



(4) 



^This value is reported with respect to the resting membrane potential. 



where 1^ is the trans-membrane current across compartment c, 
Tc is the distance of compartment c from the recording elec- 
trode, and pc is the extracellular resistivity at compartment c. 
The value of pc for soma and proximal apical dendritic compart- 
ments is 643 ^2 cm (Lopez-Aguado et al., 200 1 ) , while for the basal 
and distal apical dendritic compartments it is 260 and 287 Q cm, 
respectively (Lopez-Aguado et al., 2001). 

RESULTS 

NETWORK SIMULATIONS 

The network is simulated in the absence and in the presence of 
astrocytic activity, i.e., a feedback to neurons in respect to synap- 
tic transmission. Figure 2 shows representative action potential 
(AP) bursts in the axon compartment of the CA3 neuron that lead 
to subsequent firing in the CAl neurons 1-4, with a slight delay, 
as expected from the flow of the input signal in the CAl region. 
The network was then simulated again, while incorporating the 
astrocytic response to neuronal synaptic transition, i.e., gliotrans- 
mission, that leads to an induction of SICs in the CAl neurons. 
We show a representative astrocytic response and its associated 
effects on the firing of the CAl neurons (Figure 3). It is apparent 
by comparing the two different simulations (compare Figures 2 
and 3; and Figures 4A and 4B) that the astrocyte makes the CAl 
neurons to fire in synchrony, and in a different pattern. Note that 
the astrocyte induces SIC twice (data not shown; this effect is 
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captured by the gating variable h in Equation 3) while the astro- 
cytic Ca^+ concentration is above the threshold for triggering 
exocytotic release of glutamate (see the shaded area in Figure 3). 

LOCAL FIELD POTENTIALS 

Local field potential (LFP) calculation at the point of interest is 
based on linear source approximation in the absence and pres- 
ence of the astrocyte. It is apparent from the peaks of LFP that the 
CAl neurons fire synchronously under the astrocyte-dependent 
pathway (Figures 4 and 5). Hence, after the addition of the 
astrocyte-induced SICs, the CAl neuron pack starts firing at a 
different rate, which is dependent on astrocyte Ca^+ oscillations 
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FIGURE 2 I Model simulations of the CA3-CA1 pyramidal neuron 
network (shown in Figure 1), in the absence of s.o. astrocytic activity, 
showing action potential generation in the axon compartment of the 
CA3 and CAl pyramidal neuron models. The CA3 neuron was given a 
constant current of 0.6 nA , in the soma compartment, which evokes a 
high-frequency burst of ~100Hz and of various durations, which repeats 
itself at ~0.66 Hz. There were corresponding responses in the 
post-synaptic CAl neurons. 



(Figure 4). The LFP duration is much shorter in the absence of 
astrocytic activity (Figure 5A) compared to the LFP duration in 
the presence of an active astrocyte (Figure 5B); the input signal, 
from the CA3 neuron, was same for both pathways. To assess 
a possible difference in the dominant frequency between the 
two pathways, we used Fast Fourier Transform (FFT) spectrum 
analysis (Figures 5D,E). Based on this analysis, the astrocyte- 
independent pathway has wide-spread activity with multiple 
peaks, a pattern likely resulting due to the differences in the firing 
rate of each CAl neuron within the pack; the main proportion 
of the activity forms an envelope, around 4-1 1 Hz, lacking any 
typical rhythm (see Figure 5D). The FFT analysis of the astrocyte- 
dependent pathway shows a shift towards lower frequency with 
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FIGURE 4 I Calculated field potential based on model simulations 
of the network in the absence (A) and in the presence (B) of the 
astrocyte. The simulation omits the first 10s of simulation as, in the 
model, astrocytic Ca^+ oscillations begin after 10 s of the CA3 
neuronal activity. 
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FIGURE 3 I IVIodel simulations of the CA3-CAT pyramidal neuron 
network (shown in Figure 1), in the presence of s.o. astrocytic activity 
i.e., gliotransmission. The CAS neuron was stimulated again in the soma 
with an input current of 0.6 nA (therefore the CA3 neuron activity is not 
shown). Representative astrocytic Ca^+ cycle and its effect on the CAl 
neurons are shown. The shaded area of the astrocytic Ca^+ indicates the 
time during which SICs are induced in the basal dendrites of the CAl 
neurons 2-4 and the apical dendrites of the CAl neuron 1. 
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FIGURE 5 I Astrocytic contribution to field potential recordings. 

Duration of evoked field potentials in the absence of astrocyte (A), in the 
presence of astrocyte activity (B), and in the presence of astrocytic activity 
along with an enhanced eNMDAR conductance in CAl neurons (C). 
Corresponding FFT spectra are shown in (D,E, and F), respectively. Note 
that FFT spectra are calculated based on complete simulated data of 50 s. 
Asterisks denote the effect of astrocyte-induced SICs. 
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sharpened tuning, resulting in an occurrence of a delta rhythm 
(Figure 5E, asterisk), possibly due to synchronous firing of the 
CAl neurons in tune with the astrocyte gliotransmission. 

EPILEPTIFORM SIMULATIONS 

NMDARs are known to play a role in the initiation and propaga- 
tion of epileptiform discharges (Auzmendi et al., 2009). To study 
the possible effect of glutamatergic gliotransmission on the eNM- 
DARs, as demonstrated in biological preparations [e.g., (Fellin 
et al, 2004)], we simulated the astrocyte-dependent pathway 
after doubling the eNMDAR conductance in the CAl neurons 
(Table 1). This computational alteration emulates a biological 
condition of an increased presence of plasmalemmal eNM- 
DARs (due to an overexpression/increased trafficking) and/or 
flux through an unaltered number of channels. As a control, 
we simulated the same network under conditions incapacitating 
gliotransmission of the s.o. astrocyte, caused by clamping astro- 
cytic Ca^+ at 100 nM, which is below the threshold for glutamate 
release. The simulation data suggest that the presence of the astro- 
cyte lacking gliotransmission does not affect the firing of the 
CAl neurons overexpressing eNMDARs upon their stimulation 
by the CA3 neuron (compare Figure SI and Figure 2). However, 
the astrocytic activity/gliotransmission may induce epileptiform 
discharges in the CAl neurons via an increased eNMDAR activ- 
ity, as evidenced by the grossly exaggerated duration of LFPs 
(Figure 5C). FFT spectrum analysis of the epileptiform activity 
further indicates that the astrocyte activity tunes the CAl neurons 
to fire in a delta rhythm, while other surrounding frequencies 
appeared suppressed (Figure 5F, asterisk). 

DISCUSSION 

Herein we presented a biophysically detailed mathematical frame- 
work of the CA3/Schaffer coUateral-CAl pyramidal neuron signal 
transduction that is based on experimental reports suggesting 
a possible role of astrocytes in contextual memory retrieval by 
synchronizing CAl pyramidal neuronal network and achieving 
a delta rhythm. We were interested in studying a small neu- 
ronal network for two reasons: (i) a small network model can 
be kept biophysically detailed, and (ii) to overcome the experi- 
mental limitations in biological preparations (i.e., a lack of the 
spatial resolution to record field potentials from such a small 
network, as opposed to recording from a much larger area in 
reality). We did analysis of the simulated local field potentials gen- 
erated in the absence and in the presence of astrocytic activity. 
The results suggest that astrocytic gliotransmission potentiates 
the synapses through eNMDARs, synchronizes CAl neurons and 
modulates their overall firing rate. Frequency domain analysis of 
the simulated LFP suggests that astrocytic Ca^+ signaling shifts 
the neuronal activity in the CAl region from a broad spectrum 
frequency lacking any typical rhythm towards a delta frequency. 
Interestingly, a recent study (Grover et al, 2009) described the 
induction of long-term potentiation in the CAl region due to a 
burst stimulation over a broad frequency range centered around 
delta frequency, albeit the underlying mechanism was not inves- 
tigated. Thus, our in silico model offers a possible computational 
explanation to this delta frequency and it should be additionally 
explored in biological preparations. 



Delta rhythm is a low frequency brain activity which is impli- 
cated in cognitive domain such as attention, salience detection, 
and subliminal perception (Knyazev, 2012). Moreover, it is 
known that sensory cues innervate the CAl pyramidal neurons 
associated with the context. Under such circumstances repeti- 
tion of an item stimulus will recover the contextual state, in the 
entorhinal cortex, in response to firing of the associated pyrami- 
dal cells in the CAl region. In our circuitry model, we observe that 
astrocytic input, which makes CAl pyramidal neurons to fire in 
a synchronous manner, eliminates group firing activity at other 
frequencies except delta, suggesting a possible role of astrocytes 
in memory retrieval. From a functional point of view, this finding 
proposes a cognitive role of astrocytes where an environmental 
cue demands attention, e.g., delta rhythm generation in visual 
cortex depends on the activity to which attention is paid (Lakatos 
et al., 2008). The model suggests that similar activity of delta 
rhythm generation might exist in the CAl region of hippocampus 
during events of contextual memory retrieval. 

NR2B subunit overexpression enhances chronic pain in mice 
(Zhuo, 2009). Additionally, NR2B-containing NMDARs are 
linked to cell death and seizures (Ding et al., 2007). To test the 
functional consequence of NR2B overexpression on CAl neu- 
ronal activity, we increased the conductance of NR2B-containing 
NMDARs which is equivalent to an insertion of eNMDARs in 
the plasmalemma or enhanced unitary conductance of the eNM- 
DARs. Such a condition simulates an increased access of astrocytic 
transmission to the NR2B-containing NMDARs. The results sug- 
gest that enhanced conductance of NR2B-containing NMDARs 
leads to a high-frequency oscillation of field potentials. The dura- 
tion of these field potentials is ~600ms which is equal to the 
duration of SIC via NR2B-containing NMDARs. This activity 
could be a molecular correlate of seizures seen in the hippocam- 
pus, which are sensitive to NR2B subunit antagonists (Kolarova 
et al, 2003; Butler et al, 2010). 

In our present work, we have not adjusted our model for 
micro -anatomical restructuring at the tripartite synapse, which 
would almost inevitably occur due to the dynamics at astrocytic 
filopodia/perisynaptic processes (Hirrlinger et al., 2004; Lavialle 
et al, 2011). Such restructuring could lead to changes in the 
proximity of neuronal and astrocytic release and/or input sites. 
Future dynamic modeling to simulate micro-anatomical arrange- 
ments would be needed to computationally address the role that 
astrocytic structural plasticity may play in tuning neuronal firing 
rhythm. 
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